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ABSTRACT 

We present estimated redshifts and derived properties for a sample of 1663 galaxies 
with Ks < 22 (Vega), selected from 50.4 arcmin^ of the GOODS/CDFS field with deep 
ISAAC imaging, and make an extensive comparison of their properties with those of 
the extremely red galaxies (ERGs) selected in the same field. We study in detail the 
evolution of Kg-selected galaxies up to redshifts 2; ~ 4 and clarify the role of ERGs 
within the total Kg-band galaxy population. Wc compute the total Kg-band luminosity 
function (LF) and compare its evolution with the ERG LF. Up to {zphot) — 2.5, the 
bright end of the Kg-band LF shows no sign of decline, and is progressively well 
reproduced by the ERGs with increasing redshift. We also explore the evolution of 
massive systems present in our sample: up to 20%-25% of the population of local 
galaxies with assembled stellar mass M > 1 x IO^^Mq has been formed before redshift 
z ~ 4, and contains ~ 45% to 70% of the stellar mass density of the Universe at that 
redshift. Within our sample, the comoving number density of these massive systems 
is then essentially constant down to redshift z 1.5, by which point most of them 
have apparently evolved into (I — Ks)-selected ERGs. The remaining massive systems 
observed in the local Universe are assembled later, at redshifts z < 1.5. Our results 
therefore suggest a two-fold assembly history for massive galaxies, in which galaxy/star 
formation proceeds very efficiently in high mass haloes at very high redshift. 

Key v^rords: galaxies: evolution - galaxies: formation - galaxies: high-redshift - 
galaxies: luminosity function, mass function 



1 INTRODUCTION 

The study of near-infrared (near-IR) galaxies, in particular 
K-band selected galaxies, has increasingly been recognised 
as an efficient method for setting constraints on the forma- 
tion epoch of massive galaxies. The rest-frame K-band is 
the most suitable photometric tracer of stellar mass, rel- 
atively independent of the galaxy star formation history. 
Thus, observer-frame K-band selected samples are expected 
to contain all of the most massive galaxies at least up to red- 
shifts 2: ~ 1 — 2. Based on the analysis of near-IR-selected 
galaxy samples, different authors have studied the evolution 
of the stellar mass content of the Universe up to redshift 
z ~ 2 (Glazebrook et al. 2004, Drory et al. 2004, Fontana 
et al. 2004), and in small-area surveys up to redshift 2; ~ 3 
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(Fontana et al. 2003, Dickinson et al. 2003, Rudnick et al. 
2003). These studies have provided evidence for the exis- 
tence of massive systems up to at least these high redshifts. 

Among the K-band selected sources, extremely red 
galaxies (ERGs) are commonly considered to trace the pro- 
genitors of the local early- type population at high redshift. 
In particular, in a bright sample of ERGs, Daddi et al. (2000) 
determined the existence of strong clustering, a characteris- 
tic property of local elliptical galaxies. This strong cluster- 
ing was later confirmed to extend to fainter magnitudes by 
Roche et al. (2002) and Roche, Dunlop & Almaini (2003) 
(hereafter R03). However, it is not yet clear whether the 
ERG population contains all the progenitors of the local 
early-type galaxies, or whether a total K-band sample of 
galaxies is necessary to select all of these progenitors at 
high redshift. To investigate the evolution of K-band se- 
lected galaxies and to assess the role of ERGs within the 
total K-band galaxy population, we present here the study 
of a sample of 1663 galaxies with Ks < 22 (Vega), and com- 
pare their properties with those obtained for a sample of 
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198 ERGs selected in the same field. The selection of the 
ERG and the total Ka-band samples from the same field is 
essential to accurately compare the properties of both popu- 
lations, by minimising the potential existence of systematic 
errors due to cosmic variance. 

The Great Observatories Origins Deep Survey 
(GOODS) (Giavalisco et al. 2004) is providing unprece- 
dented multiwavelength data in ~ 320 arcmin'^ centred 
on the Chandra Deep Field South (CDFS) and Hubble 
Deep Field North (HDFN). Within the GOODS/CDFS 
field, RQ3 selected 198 ERGs with Ks < 22 (Vega) and 
(I775 — Ks) > 3.92 from 50.4 arcmin^ for which deep near-IR 
data have been obtained with the Infrared Spectrometer 
and Array Camera (ISAAC) on the 'Antu' Very Large 
Telescope (Antu-VLT) (GOODS/EIS vO.5 release). This is 
the deepest significant sample of ERGs selected to date. 
The region covered by the ISAAC vQ.5 imaging within the 
GOODS/CDFS field is shown in Figure [T] The redshift 
distribution and other derived properties of the R03 ERG 
sample have been studied by Caputi et al. (2004); hereafter 
C04. Both R03 and C04 are complementary to other 
previous studies of shallower samples of ERGs in the 
CDFS (Cimatti et al. 2002, Moustakas et al. 2004). CQ4 
constructed the ERG luminosity function (LF) and deter- 
mined that its bright end does not evolve from redshifts 
{zphot} = 2.0 to (zphot) = 2.5, connecting this fact with the 
presence of progenitors of local L > L* galaxies at redshifts 
Zphot > 2. Also, C04 showed that the comoving densities 
of ERG progenitors of local L > L* galaxies at different 
redshifts are below the total expected values, indicating 
either that ERGs cannot account for all the progenitors of 
massive galaxies at redshift 2: = 0, or that the underlying 
assumption of passive evolution is not completely valid. 

Within the GOODS program, near and mid-IR data 
from the Spitzer Space Telescope will soon become avail- 
able for the GOODS/CDFS. Spitzer is currently imaging 
this field with the Infrared Array Camera (IRAC) at 3.6/^m, 
4.5/im, 5.8/im and 8.0/^m, and with the Multiband Imaging 
Photometer (MIPS) at 24/im. The availability of these data 
will provide a unique opportunity to follow up ISAAC and 
HST sources in the CDFS, and observe directly at wave- 
lengths which map the rest-frame Kg-band at high redshift. 
Also, the existence of mid-IR images will allow the study of 
the light re-emitted by dust in optical and near-IR sources. 
The analysis of mid-IR images will be an important tool to 
understand, for example, the extent to which the presence 
of dust is responsible for the extremely red colours observed 
in ERGs. 

In the present work we use the public reduced 
GOODS/EIS vO.5 data to select a full sample of galaxies 
with Kg < 22 (Vega) in the same area analysed by RQ3. 
The aim is to study the evolution of the total Kg-galaxy 
population with redshift and to compare the results with 
those obtained for the R03 ERG sample by C04. The layout 
of this paper is as follows. First, in Section 2, we give de- 
tails on the sample selection and the multiwavelength pho- 
tometry. In Section 3 we explain the photometric redshift 
techniques applied, and discuss the calibration of the red- 
shift estimates with spectroscopic data and with the previ- 
ous estimates for ERGs determined by C04. In Section 4 we 
present our results, and compare with the results for ERGs 
obtained by C04: we discuss raw and dust-corrected Hubble 



-27:42:00 



U 
M 
Q 



-27:48:00 



-27:54:00 



_L 



53:12:00 53:06:00 53:00:00 

RA (deg.) 

Figure 1. Schematic diagram of the sub-region covered by deep 
VLT/ISAAC near-infrared imaging in the GOODS/EIS vO.5 re- 
lease (shaded area), within the larger GOODS/CDFS field with 
full HST-ACS coverage. 



diagrams, the evolution of the LF, and the evolution of mas- 
sive systems up to redshift Zphot ~ 4. In Section 5, we pre- 
dict IR magnitudes and colours for our Ks-selected galaxies 
in three Spitzer/IRAC channels: 3.6/im, 4.5^m and 5.8^m. 
Finally, in Section 6 we offer some concluding remarks. We 
adopt throughout a cosmology with Ho = 70 kms~^Mpc~^, 
^Im = 0.3 and Qa = 0.7. 



2 THE SAMPLE 

2.1 Multiwavelength photometry 

GOODS observations include optical and near-IR imaging 
in the B, V, I775 and z bands with the Advanced Camera 
for Surveys (ACS) on board the Hubble Space Telescope 
(HST), and J, H and Ks bands with the ISAAC/VLT. The 
public reduced ISAAC images (GOODS/EIS vO.5 release) 
consist of eight maps covering approximately one third of 
the whole GOODS/CDFS area with ACS/HST coverage. 
We constructed our Ks-band catalog using the reduced Ks- 
band images and extracting sources with the public code 
SEXTRACTOR (Berlin & Arnouts, 1996) on each map sep- 
arately. As one of our aims was to compare our results with 
those obtained for ERGs by R03 and C04, we needed to se- 
lect the Ks-band sample in the same way and, consequently, 
we used the same SEXTRACTOR parameters used by R03 
for the extraction of sources, i.e. a detection threshold of 
lAasky in at least six contiguous pixels. We employed a 
Gaussian kernel with a 3-pixel FWHM and the correspond- 
ing weighting maps. The low signal-to-noise edges of each 
map have been excluded from the detection procedure, leav- 
ing a total effective area for detections of 50.4 arcmin^. There 
is some overlap between the eight Ks-band maps (even after 
the trimming of the edges). Thus, as we ran SEXTRAC- 
TOR on each of the maps separately, we discarded repeated 
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sources appearing on different map catalogs. We used to- 
tal Kron-type magnitudes for limiting the sample to objects 
with magnitude Ks < 22 (and above a 3a detection limit, 
although only a few objects with Ks < 22 were below this 
threshold). The number of sources in this first Ks = 22 cat- 
alog was 1756. For measuring colours and applying photo- 
metric redshift algorithms, we used circular 2-arcsec diam- 
eter aperture magnitudes. This aperture size is the same as 
used for the ERGs by R03 and C04. We ran SEXTRAC- 
TOR in 'double-image mode' to perform aperture photom- 
etry on the J and H band images, using the corresponding 
Ka-band images for the detection of sources. For a minority 
of objects (<10%), a 2-arcsec diameter aperture produced a 
difference > 0.5 mag between the aperture and the total Ks 
magnitudes. For these objects we used a larger aperture of 
2^?5i 2.83 -arcsec diameter. The aperture sizes we adopted 
for the near-IR photometry have been deliberately selected 
to coincide with aperture sizes used to construct the ACS 
catalogs, as we explain below. 

We looked for counterparts of our 1756 Ks-band sources 
in the public r.l.Oz GOODS/CDFS ACS catalogs, within 
a matching radius of 1 arcsec. The GOODS/ ACS catalogs 
have been constructed running SEXTRACTOR on the ver- 
sion vl.O of the reduced stacked GOODS/CDFS ACS im- 
ages. The extraction of sources has been made on the z-band 
images and SEXTRACTOR has been run in 'double-image 
mode' to perform photometry on the B, V and I775 bands. 
The GOODS/ ACS catalogs provide photometry in 11 dif- 
ferent circular apertures for each source. For each of our 
Ks-band counterparts, we used the magnitudes correspond- 
ing to the same aperture size we had selected in the near-IR 
bands. We found that 25/1756 of our Ks < 22 sources did 
not have a counterpart in the GOODS / CDFS ACS catalog. 
We made an individual inspection of each of these 25 sources 
on the Ks-band images and found that 13/25 were spuri- 
ous Ks-band detections, most of them close to the (already 
trimmed) edges of the corresponding Ks-band map and, in 
some cases, with no flux in the J or H bands. Another one of 
the 25 sources appears in a region where the ACS images are 
damaged by satellite tracks, so all the optical information is 
missed for this object. After exclusion of these 14 sources, 
the remaining sample had 1742 objects. In each case of non- 
detected flux (SEXTRACTOR magnitude 99.0) in any op- 
tical or near-IR band, we measured the flux manually in the 
corresponding aperture using the IRAF task 'phot' on the 
stacked GOODS HST/ACS Vl.O and the VLT/ISAAC vO.5 
images. For faint objects in every ACS band (all of them 
ERGs), we repeated the same procedure explained in C04, 
leaving only a few sources as formally non-detected in the 
V or redder bands in the multiwavelength catalogs we used 
as input for the photometric redshift algorithms. 

As a final comment on the photometry of Ks-selected 
galaxies, we note that we used the photometric measure- 
ments made on the public GOODS Hubble Ultra Deep Field 
(HUDF) images in the B, V, I775 and z bands for three 
ERGs. These objects were initially recognised as possible 
z > 4 candidates, but the better quality of the HUDF data 
allowed to reveal them as being at lower redshifts (C04). 



2.2 Completeness 

We used the IRAF tasks 'gallist' and 'mkobjects' to cre- 
ate a Ks < 22 mock catalog of 200 objects with a power- 
law luminosity distribution and insert them into each of the 
eight ISAAC Ks-band maps. We ran SEXTRACTOR again 
on these images, using the same original extraction param- 
eters, and checked the fraction of artificial sources recov- 
ered in each case. From this, we estimated that our sam- 
ple was 100%, 89%, 86%, and 80% complete at Ks < 20.0, 
20.0 < Ks < 21.0, 21.0 < Ks < 21.5 and 21.5 < Ks < 22.0, 
respectively. 



2.3 The ERG sample 

In this work, we compare the redshift distribution and de- 
rived properties for Ks < 22-selected galaxies with the cor- 
responding results obtained by C04 for the ERGs selected 
by R03 in the same field. However, in order to assess the 
validity of the comparison of the present Kg-selected sample 
with a subsample of ERGs obtained in an independent way, 
we checked the number of objects with (I775 — Kg) > 3.92 
among the new Ks < 22 sources, obtaining a number of 205. 
186/205 of these objects are in common with the sources se- 
lected by R03 as ERGs with the same colour cutoff (which 
are 198 objects in total), producing an overlap of both ERG 
samples >90 %. Thus, the comparison of results for the 
present Ks-selected sample and the R03 ERGs is entirely 
valid. 



2.4 Star/galaxy separation 

We cleaned our catalogs for stars/QSOs using the 
SEXTRACTOR stellarity parameter 'CLASS_STAR', as 
measured on the ACS z-band images and appears 
on the GOODS/CDFS ACS catalogs. We found that 
107/1742 Ks < 22 sources had z-band counterparts with 
CLASS_STAR> 0.8. 64 out of these 107 objects are con- 
firmed by the SIMBAD astronomical database as stars or 
QSOs, so we removed them from our catalogs, leaving 1678 
sources for the photometric redshift algorithm input files. 
The inspection of the redshift output catalogs showed that 
15 out of the remaining 43 objects with CLASS_STAR> 0.8 
could not be satisfactorily modelled by any galaxy spectral 
energy distribution (SED) template (cf. Section 3) and had 
zero probability of being at any redshift. Thus, we considered 
that these 15 objects were also very likely to be stars/QSOs 
and excluded them from our sample, leaving 1663 sources in 
the final Kg < 22 galaxy catalog. 



2.5 Number counts 

Figure|5|shows the differential number counts for our sample 
of Ks-selected galaxies (empty circles) in comparison with 
the number counts for ERGs obtained by R03 (filled circles) . 
Both sets of counts have been corrected for incompleteness. 
The ERG number-count curve is characterised by a flatten- 
ing of the count slope beyond Ks « 19.5. In Section 4.3, we 
show that this fact is related to the existence of a turnover 
in the ERG LF at redshift Zphot ~ 1- 
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Figure 2. Differential number counts for our Kg-selected galaxies 
(empty circles), in comparison with the number counts for the 
ERGs selected in the same field by R03 (filled circles). 



that neither the HYPERZ primary nor secondary solutions 
provided as low a redshift estimate as the BPZ determina- 
tion. In these cases, we adopted the BPZ estimates as the 
definitive redshifts. One extra source had a very high redshift 
HYPERZ primary solution and a low redshift secondary so- 
lution in agreement with BPZ. However, the real redshift for 
this source is Zphot = 1.669, based on spectroscopic deter- 
minations (cf. Section 3.2), which both photometric redshift 
algorithms fail to correctly estimate. In this case only, we 
adopt the spectroscopic value as the definitive redshift. For 
the eight sources with either BPZ or spectroscopic redshift 
replacements, we do not have the complementary informa- 
tion provided by HYPERZ in the output, i.e. SED best-fit 
template, k-corrected Ks-band absolute magnitude, age, etc. 

As we explained in Section 2.1, after running HYPERZ 
on our catalog of 1678 Ks < 22 sources, we determined that 
15 of these objects were likely to be stars/QSOs, based on 
their stellarity parameter and the fact that they did not have 
a satisfactory fitting by any of the templates in the HYPERZ 
SED library. After rejecting these objects, our final Ks < 22 
catalog consisted of 1663 galaxies, including ~ 95 % of the 
ERGs in the R03 sample. (The few ERGs mentioned in C04 
as likely stellar contaminants are now automatically rejected 
from the present Kg < 22 galaxy sample). 



3 REDSHIFT ESTIMATES 

3.1 Photometric redshift techniques 

We used the pubhc code HYPERZ (BolzoneUa, Miralles & 
Pello, 2000) to compute photometric redshifts {zphot) for our 
sample of Kg-selected sources, using the seven passbands de- 
scribed in Section 2.1 (B, V, I775, z, J, H and Ks) and the 
GISSEL98 SED template library (Bruzual & Chariot 1993). 
A detailed explanation of the HYPERZ performance has 
been given in C04. We only note here that we refer to the 
HYPERZ 'primary solution' as the redshift estimate Zphot 
(and corresponding set of parameters) which produce the 
maximum likelihood in the SED fitting. The 'secondary so- 
lution' is the second most likely solution in parameter space. 
We ran HYPERZ using the same parameter values as in C04. 
In particular, we applied a Calzetti et al. (2000) reddening 
law with the V-band extinction {Ay) varying between and 
1. For those objects taking the maximum possible reddening, 
we made a second HYPERZ run allowing the V-band extinc- 
tion Av to vary from to 3. We used the public 'Bayesian 
photometric redshift' (BPZ) code by Benitez (2000) to ob- 
tain a second, independent redshift estimate for a subset 
of the sources that HYPERZ found to be at high redshift, 
namely: 

• very bright Ks-band sources with primary solution 

Zphot ^ 2 

• sources with significant flux in the B-band (Bab < 27.5) 
and primary solution Zphot > 3 

• all sources with primary solution Zphot > 4. 

For those objects which BPZ also placed at high redshift, 
we accepted the HYPERZ primary solution as the definitive 
solution for the object. On the other hand, when BPZ deter- 
mined a lower redshift estimate, we adopted the HYPERZ 
secondary solution, which was in agreement with BPZ in 
most of these alternative cases. For seven sources, we found 



3.2 Photometric redshift calibration 

We used the public spectroscopic redshifts available for 
the GOODS/CDFS to calibrate our photometric redshift 
estimates. Spectroscopic data in the GOODS/CDFS 
have been obtained by different authors (Cimatti et al. 
2002b, Vanzella et al. 2004, Le Fevre et al. 2004, Szokoly 
et al. 2004, among others). All the spectroscopic red- 
shifts available for sources in the r.l.Oz GOODS/CDFS 
ACS catalog have been compiled in a master list at 
http: / / www.eso.org/science /goods /spectroscopy / 
CDFS JVIastercat. The spectroscopic redshifts ob- 
tained as part of the K20 survey (Cimatti 
et al. 2002b) have recently made public at 
|http:/ /www.arcetri.astro.it/~k20 /spe_release_dec04/index| 
htmL 

Figure EI a) compares our redshift estimates with the 
spectroscopic redshifts for the Ks < 22 sources included in 
different CDFS spectroscopy samples. We restrict the com- 
parison to sources with good quality flags for the spectro- 
scopic redshift determinations. The error bars in the photo- 
metric redshifts correspond to Icr-confidence levels. We ob- 
serve quite good agreement between the redshift estimates 
and the real redshifts in most cases. However, there is an 
excess of low redshift sources which we identified as being 
at redshifts Zphot > 2.5. These sources with mismatching 
results are indicating that the photometric redshift algo- 
rithms might have a tendency to favour high redshift val- 
ues when degenerate solutions in parameter space exist. To 
compensate for this effect, we calibrate our photometric red- 
shifts with the following criteria: for sources with signifi- 
cant B-band flux (Bab < 27.5) and a large difference be- 
tween the HYPERZ primary and secondary solutions, i.e. 
{zprim — Zsec) > 1, we adopted the HYPERZ secondary 
lower-redshift solution in the cases in which this redshift 
had at least 40% of the probability of the primary one. For 
sources with Bab > 27.5 and {zprim — Zsec) > 1, which are 
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Figure 3. Photometric redshifts vs. spectroscopic redshifts for 
Ks<22-selected galaxies in the GOODS/CDFS deep ISAAC field: 
a) before the calibration; b) after the calibration. 



more likely to be at high redshift, we adopted a threshold 
of 70% between the relative probabilities of the secondary 
and the primary solutions to adopt the former as the pre- 
ferred answer. These criteria have been applied in addition 
to those explained in Section 3.1, without affecting the cases 
in which the HYPERZ secondary solution was already se- 
lected for providing a better agreement with the correspond- 
ing BPZ estimation. In summary, we adopted the HYPERZ 
secondary solution for 136/1663 sources, while 1519/1663 
remained with the primary one (and 8/1663 had either BPZ 
or spectroscopic value replacements) . The comparison of our 
photometric estimates with the spectroscopic redshifts after 
the calibration procedure is shown in Figure l^b). The me- 
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dian of the relative errors between our estimates and the real 
redshifts is \zphot^Zspc.c\/{i + Zspec) ~ 0.05. In particular, we 
obtain a median of \zphot — Zspec\/{i + Zspec) = 0.07 for the 
compared sources with Zphot > 1.5. These values are among 
the typical uncertainties obtained in photometric redshift 
estimates available in the literature. 

Optical U-band data in the CDFS has been obtained by 
Arnouts et al.(2001). Given the shallower depth and poorer 
resolution (1.5") of the U-band images in comparison with 
those corresponding to the GOODS/HST-ACS, we decided 
not to include the former in the input for the photometric 
redshift algorithms. However, for consistency, we checked for 
the presence of counterparts of our high-redshift Ks-selected 
galaxies among the CDFS U-band sources. We found that, 
as expected, none of our Ks-selected galaxies with estimated 
redshift Zphot > 2.8 has a counterpart in the CDFS U-band 
source catalog within a matching radius of 1 arcsec. Nev- 
ertheless, it should be kept in mind that fainter U-band 
counterparts could exist in some cases. Thus, the lack of 
deep U-band data could potentially produce that a fraction 
of our high-redshift {zphot > 2.8) sample is contaminated 
with faint low-redshift galaxies with degenerate solutions in 
photometric redshift space. 

On the other hand, while we were performing the 
present study, a new release (vl.O) of the GOODS/CDFS 
ISAAC data in the Ks and J bands has been made public. 
The new release includes the determination of a uniform 
photometric calibration across the entire GOODS/CDFS 
field. We revised the Ks and J-band photometry of our Ks- 
selected sample in the new ISAAC maps. After running HY- 
PERZ and BPZ under the same conditions as explained in 
Section 3.1, we determined the existence of a basically neg- 
ligible impact of any possible difference between the ISAAC 
vO.5 and vl.O calibrations on our estimated redshifts. 



3.3 The ERG photometric redshifts 

The calibration procedure explained in Section 3.2 was not 
implemented in C04 because CDFS spectroscopic redshifts 
were not yet completely available and because the overlap 
with the R03 ERG sample was small. Thus, we need to ver- 
ify that the present selection of a Ks-band sample and the 
calibration of the photometric estimates based on spectro- 
scopic redshifts still produces results in agreement with C04 
for the ERG subsample. 

Figure |1] shows the new redshift estimates we deter- 
mined for the ERGs versus the redshift estimates obtained 
by C04. The error bars correspond to Icr-confidence lev- 
els. Both sets of redshifts appear to be in very good agree- 
ment, except for a few outliers which correspond to cases of 
very strong degeneracies in parameter space. This correla- 
tion confirms that, in spite of the differences in the redshift 
calibration implemented in the present work, the redshift es- 
timates obtained for the ERGs are consistent and the com- 
parison of derived properties for our total Ks-selected sample 
with the results obtained in C04 for the ERGs is valid. 
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Figure 4. Rcdshift estimates determined in the present worlc for 
tlie R03 ERGs versus the redshift estimates obtained in C04. 



4 RESULTS 

4.1 The redshift distribution 

Figure |SJ a) shows the redshift distribution of the total sam- 
ple of Ks-selected galaxies (solid lines) compared with the 
redshift distribution for the R03 ERGs obtained by C04 
(dashed lines). The corresponding histograms have been 
constructed taking into account the probability density dis- 
tribution in redshift space for each object, as given by HY- 
PERZ in the output. The consideration of probability den- 
sities yields more realistic representations of the redshift 
distributions, because they take into account all the un- 
certainties and degeneracies inherent to photometric red- 
shift techniques. For the minority (<9%) of objects for 
which we adopted the HYPERZ secondary solution or a 
BPZ/spectroscopic replacement, we only took into account 
a single redshift determination, as any information on a well- 
determined probability density distribution in redshift space 
was missing for these objects. Both histograms in Figure|5Ja) 
include corrections for incompleteness, i.e. we multiplied by 
a factor 1.12, 1.16 and 1.25 the contribution of sources with 
20.0 < Ks < 21.0, 21.0 < Ks < 21.5 and 21.5 < < 22.0, 
respectively. The distribution of the total Ks-selected galaxy 
sample peaks at redshift Zphot ~ 0.5 — 1.0, while the ERG 
distribution has a maximum at Zphot ~ 1.5. The differ- 
ence in the redshift distributions is a consequence of the 
ERG colour cutoff. Figure |SJb) shows the relative fraction 
of ERGs among Ks-selected galaxies, allowing for an eas- 
ier comparison between both populations. We see that the 
fraction of Ks-selected galaxies which are (I775 — Ks) > 3.92 
ERGs is approximately constant (~ 20%-30%) between red- 
shifts 1.5 < Zphot < 4.0. We only find two sources with es- 
timated redshifts Zphot > 4 in our present sample and these 
are the same two ERGs reported as confident z > 4 candi- 
dates in C04. 
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Figure 5. a) The redshift distribution of the total Ks < 22 
sample of galaxies (solid lines) compared with the rcdshift distri- 
bution of the R03 ERGs (dashed lines). The adoption of a prob- 
ability density distribution in redshift space for each source and 
incompleteness correction factors produce a non-integer number 
of sources in each redshift bin. b) The percentage of ERGs found 
among Ks-selected galaxies as a function of redshift. 



4.2 The Hubble diagram 

In Figure |H| we plot the observed Kg magnitude versus pho- 
tometric redshift Zphot for the R03 ERGs (filled circles), as 
they were obtained by C04, and for all the other Ks<22- 
selected galaxies (empty circles). The solid line shows the 
empirical K-z relation for radio galaxies obtained by Willott 
et al. (2003), which approximately corresponds to the pas- 
sive evolution of a 3L* starburst formed at redshift 2; = 10 
and indicates the behaviour of the most massive galaxies 
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Figure 6. Total raw Kg magnitudes versus pliotornetric red- 
sliifts Zphot- Tlie empty circles correspond to the Ks<22-selectcd 
galaxies which are not in the R03 ERG sample. The filled circles 
correspond to the ERGs in the R03 sample, with the redshifts 
estimated by C04. The solid line shows the empirical K-z rela- 
tion for radio galaxies as obtained by Willott et al. (2003), which 
approximately coincides with the passive evolution of a 3L* star- 
burst formed at redshift z = 10. The dotted line is a nominal 
extrapolation of the same law. 



formed at very high redshifts. The redshift associated with 
each Ks-selected galaxy in the present sample is either a 
HYPERZ solution or a BPZ/spectroscopic replacement, in 
accord vsrith the criteria explained in Sections 3.1 and 3.2. 

The comparison of the Y^s-Zphot relation for extremely 
red and bluer Kg-selected galEixies yields two main conclu- 
sions. First, ERGs show the same large dispersion in the 
Hubble diagram as all the other Ks-selected galaxies (al- 
though the ERG colour cutoff produces the absence of these 
objects at the lower redshifts). Second, as was discussed in 
C04, the Hubble diagram is characterised by a lack of ob- 
jects near the radio galaxy line at redshifts Zphot > 3. We 
see that this effect is still present when all the other Ks- 
selected galaxies arc included. However, we note the pres- 
ence of two relatively bright objects at Zphot > 3 which 
are not (I775 — Kg) > 3.92 ERGs. These bright objects at 
Zphot > 3 have (J — Kg) > 2.2 (Vega), but some of them are 
missed when a typical ERG (I775 — Kg) colour cut is applied 
(in the present case, in particular, these two sources have 
(I775 — Kg) < 3.7). These objects could correspond to galax- 
ies in the final stages of a continuous but decaying process 
of star formation, sufficiently old (~ 1 — 2 Gy) to have de- 
veloped a 4000A-break which produces red observer-frame 
(J — Kg) colours. Yet, in some of them, the still rather active 
star formation produces considerable amounts of ultraviolet 
(UV) flux which, in conjunction with only modest dust ex- 
tinction, prevents these sources from displaying extremely 
red (I775 — Kg) colours. In Section 4.4, we present some fur- 
ther discussion on the influence of the selection effects on 



the inclusion of the most massive galaxies when construct- 
ing Ks-band galaxy samples. 

4.3 The evolution of the Kg-band/ERG luminosity 
functions 

In this section, we study the rest-frame Kg-band LF for 
the ERGs and the total sample of Kg-selected galaxies at 
different redshifts, in order to investigate any possible dif- 
ference in their evolution. We extrapolated the rest-frame 
Kg-band magnitudes using the best-fit SED of each galaxy. 
These extrapolations can be large at high redshifts, but the 
k-corrections in the Kg-band are less dramatic than those 
corresponding to other bands (Poggianti 1997). Although 
the computation of a rest-frame optical LF would require 
smaller extrapolations at high redshifts, a proper compari- 
son with the ERG LF is difficult to perform at rest-frame 
optical wavelengths. The values obtained for the ERG LF 
would strongly depend on the extinction corrections Ay and, 
thus, for our purpose, a rest-frame optical LF would be more 
model-dependent than the LF in the rest-frame Kg-band. 

We computed the Kg-band LF binning the Kg < 22- 
samplc of 1663 galaxies in both redshift and absolute mag- 
nitude space, in the same way as was done for the ERG sam- 
ple by C04. A single redshift estimate has been considered 
for each source (either the HYPERZ primary, the secondary, 
or a BPZ/spectroscopic solution, as detailed in Section 3.2). 
In the cases in which we adopted the HYPERZ primary 
solution, the corresponding k-corrected absolute magnitude 
was directly obtained from the HYPERZ output. In 
all the other cases, where the estimated redshifts are rela- 
tively low, we computed the absolute magnitude Mk^ as- 
suming that the k-corrections have negligible dependence 
on the SED shape (an assumption particularly valid for the 
Kg band at relatively low redshifts). To take into account 
the limits of the survey (Kg = 22), we weighted each source 
by a factor Vmaxbin/Vmaxob3, where VmaxUn is the volume 
determined by the maximum redshift of the corresponding 
bin and Vmaxobs is the volume determined by the maximum 
redshift at which the source would still be detected in the 
survey (provided it is lower than the maximum redshift of 
the bin). We also applied weighting factors of 1.12, 1.16 and 
1.25 for sources with 20.0 < Kg < 21.0, 21.0 < Kg < 21.5 
and 21.5 < Kg < 22.0, respectively, to correct for incom- 
pleteness of the sample (Section 2.2). 
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Figure 7. The rest-frame Kg-band luminosity function at red- 
shifts (zphot) = 1-0,1.5,2.0 and 2.5 (circles, squares, up triangles 
and diamonds, respectively). The values shown are strictly above 
the completeness limits of each redshift bin. The solid line corre- 
sponds to the Schechter function fitted to the local K-band LF 
measured on 2MASS data by Kochanek et al. (2001). The dashed 
line shows the evolution of the LF at redshift z = 1, as estimated 
by Drory et al. (2003). 




Figure 8. The raw rest-frame Ks-band luminosity function at 
redshifts (zphoi) = 1.0, 1.5, 2.0 and 2.5 (no Vmax corrections, no 
corrections for incompleteness). The symbols are the same as in 
Figure ItI 
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Figure 9. The comparison of the rest-frame Kg-band LF ob- 
tained from the present study (empty symbols) with the LF ob- 
tained by Pozzotti ot al. (2003) using K20 survey data (asterisks): 
a) (zphot) = 1-0; b) (z^hot) = 1-5. 



Figure Q shows the total Kg-band LF at redshifts 
(zphot) = 1.0, 1.5, 2.0 and 2.5 (circles, squares, up triangles 
and diamonds, respectively). The values shown are strictly 
above the completeness limits of each redshift bin. The error 
bars correspond to the maximum of Poissonian errors and 
the errors due to cosmic variance, which we considered on 
average as 40% and 30% of the number counts at redshifts 
z < 2 and 2 < z < 3, respectively (Somerville et al. 2004). 
Although not all the Kg-selected galaxies are expected to 
display the strong clustering observed in ERGs (R03), we 
adopted the same cosmic variance error bars as in C04 for 
conservative reasons. We do not present the LF at redshift 
{zphot} ~ 0.5 because our surveyed volume is small at those 
redshifts and the LF is quite affected by the presence of large 
scale structure at « = 0.67 — 0.73 (Le Fevre et al. 2004; 
Vanzella et al. 2004). We show, for comparison, the local 
K-band LF obtained by Kochanek et al. (2001) using Two 
Micron All Sky Survey (2MASS) data, and its evolution to 
redshift 2: = 1 as estimated by Drory et al. (2003) (solid and 
dashed lines, respectively). Drory et al. (2003) found that 
the density of bright objects increases from redshift z = to 
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Figure 10. Compared evolution of ttie total Kg-band (empty symbols) and the ERG (filled symbols) LPs with redshift. The values 
shown are strictly above the completeness limits of each redshift bin. The error bars for each point are the maximum of Poissonian errors 
and the errors due to cosmic variance. 



z = 1. A similar trend was found by Pozzetti et al. (2003) 
from the analysis of the near-IR LF up to redshift (z) = 1.5, 
obtained using K20 survey data. Our results confirm this, 
and show no evidence for the bright end to decrease again 
up to at least redshifts (zphot) = 2.5. The same effect was 
observed for the ERG LF by C04. Figure confirms that the 
behaviour of a constant (or even possibly increasing) evolu- 
tion for the bright end of the LF is a property of the total 
Ks-selected galaxy population. At this point, we should re- 
mind the potential effects of the lack of deep U-band data 
on our estimated redshifts. As we discussed in Section 3.2, a 



non-negligible fraction of the high-redshift sample analysed 
here could potentially be contaminated with lower redshift 
sources with degenerate solutions in photometric redshift 
space. This fact could partially bias the trend we observe 
in the bright end of the LFs. At intermediate magnitudes, 
on the other hand, our results suggest a non-negligible de- 
crease in the density of objects, although wider and deeper 
surveys are necessary to better quantify this evolution. To 
demonstrate that none of our conclusions depend on the 
Vmax corrections or the corrections for incompleteness, we 
show in Figure |H1 the 'raw' Ks-band LF, i.e. the LF con- 
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Figure 11. The relative errors of the redshift estimates for the sources considered in the computation of the Ks-band LF at different 
redshifts. The redshift bins considered in a), b) and c) have width Azp^j^t = 0.5, while the redshift bin considered in d) has width 



structed without applying any of these corrections. We note 
that the analysis of the raw LF allows us to extract exactly 
the same conclusions as in the case of the corrected version 
presented in Figure[7| The comparison of our rest-frame Ks- 
band LF with that obtained for the K20 survey is shown in 
Figure ^ 

The four panels in Figure [THl compare the evolution of 
the total (empty symbols) and the ERG (filled symbols) Ks- 
band LFs with redshift. The latter is the same as presented 
in C04, although here we only show our LF values strictly 



above the completeness limits of each redshift bin. For each 
LF individually, we adopted the maximum of the Poisson 
errors and the errors due to cosmic variance. Any possible 
error produced by uncertainties in our redshift estimations is 
expected to be contained within the Poisson/cosmic variance 
error bars. Figure [TTI shows the relative error dZphot/{i + 
Zphot) in the redshift estimation of each one of the sources 
making the Ks-band LF at different redshifts. The absolute 
error dzphot corresponds to la-confidence levels, as obtained 
from HYPERZ. The sources with dzphot ~ correspond to 
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cases in which the HYPERZ secondary solution has been 
adopted for the redshift estimate. In every redshift bin, the 
median of the relative errors for the redshift estimates is 
rfZphot/(l + Zphot) < 0.06 and the number of sources with 
dZphot/{l + Zphot) > 0.2 is very small (<3% in all cases). 

A clear feature in the evolution of the LF is the increas- 
ing role of ERGs in reproducing the bright end of the total 
Ks-band LF with increasing redshift, reaching a maximum at 
{zphot) ~ 2.0, where a considerable fraction of the brightest 
Ks-selected galaxies are ERGs. Also, from the comparison of 
the total and the ERG Ks-band LFs at redshift (zphot) = 1.0, 
we conclude that the latter has a bell-like shape, consistent 
with a Schechter function with a flat slope. This result is 
in agreement with the different shapes determined for the 
LFs of red and blue galaxies at different redshifts (e.g. Lilly 
et al. 1995; Pozzetti et al 2003; Kodama et al. 2004). At 
redshift (zphot) = 1.5, the total Kg-band and the ERG LFs 
are more similar, but a significant density of blue galaxies is 
still necessary to account for the total Ks-band galaxy pop- 
ulation at intermediate magnitudes. At 2 < Zphot < 3, our 
completeness limits only allow us to explore the bright end 
of the Ka-band LF, which appears mainly dominated by the 
ERGs. 

4.4 The evolution of massive galaxies 

In C04, it was shown that the comoving densities of ERG 
progenitors of the local L > L* galaxy population, under 
passive evolution, were below the values expected by the 
evolution of dark matter haloes in a A-Cold Dark Matter 
(ACDM) formalism (Kauffmann & Chariot 1998; Somerville 
& Primack 1999). This fact has been interpreted as either 
a deficiency of the ERG population to account for all the 
progenitors of the local L > L* galaxies, or as due to the in- 
complete picture described by a passive evolution scenario. 
Our aim here is to further test these two possibilities using 
a similar approach: studying the evolution of the comoving 
number density of massive systems and, in general, the evo- 
lution of the stellar mass content of the Universe with red- 
shift. The evolution of the stellar mass density has recently 
been studied in different fields up to redshift 2 ~ 2 (Glaze- 
brook et al. 2004, Drory et al. 2004, Fontana et al. 2004), 
and in small-area surveys up to redshift z ~ 3 (Dickinson et 
al. 2003, Fontana et al. 2003, Rudnick et al. 2003). Here we 
extend the analysis of the comoving mass density evolution 
up to redshift z ~ 4, as obtained from the study of our deep 
sample of Ks-selected galaxies in the GOODS/CDFS field. 

We computed the rest-frame Kg-band luminosity of each 
of our galaxies using the k-corrected absolute magnitude 
Mk, , associated with the HYPERZ primary or secondary 
solution, on a case-by-case basis. In the cases in which we 
adopted the HYPERZ secondary solution, the absolute mag- 
nitude has been extrapolated as we explained in Section 
4.3. To estimate the mass of our galaxies, we used, in each 
case, the galaxy luminosity and a mass-to- light (M/Lks) ra- 
tio depending on the corresponding best-fit HYPERZ SED 
type and age. For the eight objects with BPZ/spectroscopic 
redshifts, any information on the galaxy SED type or age is 
missing, so these sources are not taken into account in all 
the analysis made in this section. We used the public code 
GALAXEV (Bruzual & Chariot 2003) to construct a grid of 
mass-to-light (M/Lks) ratios depending on galaxy age. In 



the rest-frame Ks-band, these ratios have little dependence 
on either the dust corrections or the galaxy star formation 
history. Thus, to model the M/Lks ratios as a function of 
age, we used only two galaxy templates, sufRciently rep- 
resentative of the different HYPERZ library SEDs. Both 
GALAXEV templates had a solar metallicity and no dust 
corrections, but corresponded to two different exponentially- 
declining star formation histories: fast (r ~ 0.1 Gy) and 
slow (r ~ 5Gy). The final M/Lks value adopted for each 
galaxy was the one corresponding to the GALAXEV tem- 
plate most similar to the best-fit HYPERZ SED, and to the 
best-fit age, obtained by interpolating between the M/Lk^ 
values in our age grid. The sufficient sampling of our age 
grid makes that the interpolation in age does not constitute 
a significant source of uncertainty. The interpolation in star 
formation history, on the contrary, might constitute a non- 
negligible source of uncertainty on the estimated masses. 
However, for the two adopted star formation histories, the 
M/Lkb ratios only differ by a factor < 2 at any given age. 
These differences are considerably smaller than those ob- 
tained with optical M/L ratios, even with the use of more 
refined grids. On the other hand, it should be noted that 
the uncertainities introduced in our mass estimates by the 
use of extrapolated luminosities in the rest-frame Ks-band 
are not larger than the uncertainties produced linking mass 
and light in rest-frame optical bands: in both cases, the pre- 
cision in the estimated mass for each galaxy relies on the 
precision of its SED modelling. We adopted a Salpeter ini- 
tial mass function (IMF), for ease of comparison with most 
of the values quoted in the literature. To compute comov- 
ing densities, we used the complete probability distribution 
in redshift space to determine the potential contribution of 
each object at different redshifts (note, however, that a more 
rigourous procedure should take into account a probability 
density in the whole parameter space). For the sources for 
which we adopted the HYPERZ secondary solution, a sin- 
gle redshift has been considered, as in these cases we cannot 
recover a well-determined probability density function. 

Figure 1121 shows the comoving number densities of 
galaxies which have assembled a stellar mass M > 5 x 
lO^^Mo (a) and M > 1 X 10"Mq (b), as a function of 
redshift. The empty circles correspond to all the Ks-band 
selected galaxies, while the filled circles indicate the con- 
tribution of the R03 ERGs selected in the same field. The 
horizontal error bars indicate the binning in redshift space. 
The vertical error bars correspond to cosmic variance (40%, 
30% and 15% of the number counts at redshifts Zphot < 2, 
2 < Zphot < 3 and 3 < Zphot < 4, respectively), which is 
the dominant source of error in this case. Vmaxbin /Vmaxoba 
and incompleteness correction factors have been applied 
throughout. In both panels (a) and (b), the star-like symbol 
(diamond) indicates the local number densities of all (early- 
type) galaxies with stellar masses above the corresponding 
mass cuts. The local number densities have been obtained 
by integrating the K-band-derived mass functions computed 
by Bell et al. (2003), converted to a Salpeter IMF. 

From inspection of Figure 1121 we see that a consid- 
erable fraction (~ 20%-25%) of the galax;ies with mass 
M > 5 X 10^°Mq and M > 1 X lO^M© is in place be- 
fore redshift Zphot ~ 4. The number density is virtually 
constant down to redshift Zphot ~ 2 and the assembly of 
the remaining massive systems appears to have been pro- 
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Figure 12. The comoving number density of galaxies with steUar 
mass M > 5 X IO^^Mq (a) and M > 1 X lO^M© (b). The empty 
circles indicate the densities of all the Kg-band selected galaxies, 
while the filled circles show the contribution of the ERGs selected 
in the same field. In both plots (a) and (b), the star-like symbol 
and diamond indicate the corresponding local values for all and 
early-type galaxies, respectively (Bell et al. 2003). The dotted- 
dashed lines indicate the evolution of the number density of haloes 
with total mass M > 8 X lO^M© (a) and M > 2 X lO^^M© (b), 
as obtained from ACDM models of structure formation. 



duced at redshifts Zphot < 1-5. It should be noted that the 
fraction of massive systems estimated to be at z > 2 would 
be somewhat overestimated if lower redshift contaminants 
were present in the high-redshift sample. Thus, the percent- 
ages quoted here should be regarded as upper limits on the 
fraction of massive systems present at high redshifts. The 
evolution observed in the number densities of sources with 



mass M > 5 X lO^'M© and M > 1 x lO^M© suggests the 
existence of a two-fold mechanism for the construction of 
massive galaxies: while most of the systems are assembled 
at relatively low redshifts, as expected in hierarchical mod- 
els of galaxy formation, a substantial fraction is assembled 
very efficiently at very early times. 

As in C04, we compare the comoving number densi- 
ties of massive galaxies with the corresponding densities of 
dark matter haloes massive enough to host these systems, as 
they are predicted by ACDM models of structure formation. 
In Figure [T^ the dotted-dashed hues represent the expected 
evolution of the number densities of haloes with total masses 
M > 8 X 1O"M0 (a) and M > 2 X IO^^Mq (b). The halo 
mass thresholds have been deliberately selected to coincide 
with the local number densities of galaxies with stellar mass 
M > 5 X 10^°Mq and M > 1 x 10"Mq, respectively. In the 
case of a single halo occupation number, we would expect 
the evolution of the number densities of galaxies to follow 
the evolution of the corresponding number density of haloes. 
This is what we observe up to redshift Zphot ^ 1. However, 
between redshifts Zphot ~ 1 and Zphot ~ 3, the predicted 
number densities of host haloes are significantly larger than 
the observed densities of already assembled massive galax- 
ies. The latter are approximately constant up to redshift 
Zphot ~ 4, while the number of haloes decrease with redshift, 
in such a way that the number of galaxies and haloes coin- 
cide again at Zphot > 3. This comparison suggests that, at 
very high redshift, the timescale in which halo merging and 
galaxy collapse occur are quite similar and, thus, the num- 
ber density of galaxies traces the number density of haloes. 
At Zphot < 3, on the contrary, halo merging could be faster 
than galaxy collapse, explaining why the observed density 
of massive galaxies is only a fraction of the density of the 
corresponding haloes. By Zphot < 1, the assembly of massive 
galaxies is mostly complete, and the number of galaxies and 
haloes coincides again. We note that, even while we have 
an excess of galaxies with mass M > 1 x 10^^ M© with re- 
spect to the corresponding number of haloes at Zphot > 3, 
this excess is not significant and could be contained within 
the error bars if the error of 15% in the number counts were 
somewhat underestimating the cosmic variance at those red- 
shifts. Future deeper surveys will be able to test whether the 
predicted number densities of haloes is still compatible with 
the observational data beyond redshift Zphot ~ 4, and set 
tighter constraints on the high formation redshifts of the 
oldest massive systems. 

As we stated at the beginning of this section, one of our 
main aims is to assess the role of ERGs to account for the 
progenitors of the most massive systems and the validity of 
the passive evolution assumption. We make clear that, of 
course, only a subset of ERGs can potentially be the pro- 
genitors of the most massive local galaxies. The whole of the 
ERG population spans a wide variety of redshifts and masses 
(C04), and many ERGs are completely irrelevant for the 
present discussion. The comparison of the number densities 
for the massive ERGs and all the massive Ks-selected galax- 
ies in Figure a) and ll2r b'l illustrates the fact that ERGs 
trace the high mass end of the Ks-selected galaxies. How- 
ever, we observe that ERGs cannot account for all the sys- 
tems with mass M > 1 x IO'^^Mq . At redshifts Zphot ~ 1 — 2, 
ERGs constitute ~ 50%-70% of all the Ks-selected galaxies 
which have assembled a stellar mass M > 1 x lO'^^ M0 . Figure 
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I12f b') also suggests the existence of an evolutionary sequence 
between massive Ks-selected galaxies at very high redshift 
and massive ERGs. We note that the comoving number den- 
sities of Ka-selected galaxies at (^phot) =3.5 is only slightly 
larger than the number density of ERGs with similar mass 
at redshift Zp^ot ~ 1, indicating that most of the massive 
systems at very high redshift might have evolved into ERGs 
by redshift Zphot ~ 1- The extremely red colours could be the 
consequence of the ageing of the stellar populations and this 
would be evidence for passive evolution since very high red- 
shift. However, as stated in C04, the colours of some of the 
reddest ERGs can only be explained by the superposition of 
an evolved stellar population and dust. This would indicate 
that, in some of these galaxies, the passive ageing of the 
stellar populations could be interrupted by the production 
of new star formation. But, except for these possible tem- 
porary periods of additional star formation, ERGs seem to 
account for the fraction of massive galaxies at high redshift 
which have plausibly been evolving under passive evolution. 
In any case, we see that the massive galaxies which could 
have evolved passively since very high redshifts are only a 
fraction of the local value and, thus, a passive evolution sce- 
nario alone is not able to explain the construction of all the 
massive systems. 

At Zphot > 2, on the contrary, we find a significant frac- 
tion of Ka-selected galaxies with mass M > 1 X 1O"M0 
which are not (I775 — Ka)-selected ERGs. Regarding the ef- 
fects of completeness, we estimate that our Ka-selected sam- 
ple is complete for galaxies with mass M > 1 x 10^^ M© up to 
redshift Zphot ~ 4, based on the median of the k-corrections 
and the maximum possible M/Lk^ ratio at that redshift. 
The median of the k-corrections is somewhat larger when 
only the ERGs are considered, and we estimate that the 
ERG sample is complete to a mass M = 1 X 1O"M0 up 
to redshift Zph^t ~ 3. Thus, incompleteness only affects the 
comparison of the ERG and all the Ks-selected galaxy densi- 
ties in the highest redshift bin. Most of our massive galaxies 
at Zphot > 2 have (J — Ks) > 2, indicating that a red (J — Ks) 
colour cut provides a more efficient method to select massive 
systems at very high redshift than an (I775 — Ks) cut (Franx 
et al. 2003). In Section 5, we predict the IR colours of our 
Ks-selected sources in different Spitzer/IRAC channels and 
present some further discussion on the efficiency of different 
pure IR colours to select high redshift galaxies. 

Another point of interest is the analysis of how the total 
stellar mass budget is distributed at different redshifts. Fig- 
ure ^1 shows the evolution of the total stellar mass density 
(circles) with redshift, as derived from our sample of Ks- 
selected galaxies. The large and small squares indicate the 
contribution of galaxies with massses M > 1 X 10" M© and 
7 X 1O^"M0 < M < 1 X 10"Mo, respectively. (Our Ks =22- 
limited sample is estimated to be complete for galaxies with 
masses M > 1 x 10" and M > 7x lO^^M© up to redshifts 
Zphot ~ 4 and Zphot ~ 3, respectively). The symbols with a 
star at 2: = show the corresponding local values obtained 
from the integration of the local galaxy stellar mass function 
(Cole et al. 2001, Bell et al. 2003). The solid line is the evolu- 
tion of the stellar mass density as expected from the integra- 
tion of the star-formation rates derived from Sloan Digital 
Sky Survey (SDSS) data (Heavens et al. 2004). From Fig- 
ure^] we see that most of the stellar mass is assembled at 
redshifts 0.5 < Zphot < 1.5, in agreement with the tendency 
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Figure 13. The distribution of the stellar mass density as a func- 
tion of redshift, as derived from our deep sample of Ks-selected 
galaxies. The circles indicate the total stellar mass density. The 
large (small) squares are the contributions of galaxies with masses 
M > lxlO"M0 (7xlOi''M0 < M < IxlO^M©), for which our 
Ka = 22 -limited sample is estimated to be complete up to red- 
shift Zphot ~ 4 {zphot ~ 3). The symbols with a star at Zphot = 
show the corresponding local values obtained from the integration 
of the local galaxy stellar mass function (Cole et al. 2001, Bell et 
al. 2003). The solid line is the evolution of the stellar mass den- 
sity as obtained from the integration of the star-formation rates 
derived from Sloan Digital Sky Survey (SDSS) data (Heavens et 
al. 2004). 

previously observed by other authors. At (zphot) ~ 1.75, we 
find that the total stellar mass density is already ~ 25% - 
30% of the local value, a percentage somewhat higher than 
the value of ~ 13% obtained by Rudnick et al.(2003) at 
z ~ 2, using a small-area survey in the Hubble Deep Field 
South (HDFS), and more similar to the value of ~ 20% to ~ 
35% determined by Fontana et al. (2004) from the K20 sur- 
vey. Our observed value of the total stellar mass density at 
(zphot) = 1.75 is (1.37 ± 0.55) x lO^MoMpc"^ (taking into 
account cosmic variance), which is in excellent agreement 
with the value (1.45tQ;g2) x lO^MoMpc"^, as predicted by 
Fontana et al. (2004) by the extrapolation of the Schechter 
fits to the mass function of the K20 survey at these red- 
shifts. It is worthwhile to note, however, that the observed 
mass density at (zphot) = 1.75 in a survey limited at Ka = 20 
is only ~ 50% of the total value (Fontana et al. 2004), while 
a survey extended up to Ks = 22 is able to recover most 
of the mass at that redshift, as we see from the comparison 
with the assembled stellar mass obtained from the integra- 
tion of the star-formation rates derived from SDSS data. 
Beyond Zphot ~ 2, we observe only a slow decrease of the 
total stellar mass density, in spite of the increasing incom- 
pleteness for the intermediate and lower mass galaxies. We 
estimate that ~ 24% and ~ 18% of the local stellar mass 
was assembled at redshifts (zphot) ~ 2.5 and (zphot) ~ 3.5, 
respectively. A detailed comparison of our total stellar mass 
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Figure 14. The comparison of the total stellar mass densities ob- 
tained from the present study (empty circles, with error bars) with 
values previously obtained by other authors: Dickinson et al. 2003 
(down-triangles); Fontana et al. 2003 (asterisks); Rudnick et al. 
2003 (diamonds); Glazebrook et al. 2004 (up-triangles); Fontana 
et al. 2004 (squares and double squares for the observed and ex- 
trapolated values in the K20 survey, respectively) . The circle with 
a star at Zpjiot = show the local value of the total stellar mass 
density obtained from the integration of the local galaxy stellar 
mass function (Cole et al. 2001, Bell et al. 2003). The solid line 
is the evolution of the stellar mass density as obtained from the 
integration of the star-formation rates derived from Sloan Digi- 
tal Sky Survey (SDSS) data (Heavens et al. 2004). Dickinson et 
al. values correspond to the average of the values obtained with 
the 1-component and 2-component models with solar metallic- 
ity (see Table 3 in Dickinson et al. 2003). Rudnick et al. values 
are computed only on galaxies with rest-frame V-band luminosity 
Lv > 1.4 X IOIOL0 and, thus, miss a significant fraction of the 
mass. Glazebrook et al. values only correspond to galaxies with es- 
timated mass log(M) > 10.46 and are at least partially incomplete 
above redshift z ~ 1.2. It should be noted that all the previous de- 
terminations of stellar mass densities at redshifts z > 2 have been 
made on surveys with coverage areas ~ 10 times smaller than the 
area analysed in this study. All the values in this figure assume 
a Salpeter IMF and a cosmology with Ho = TOkms^'^Mpc"^, 
Om = 0.3 and Ca = 0.7. 



densities at different redshifts with other values found in 
the literature is presented in Figure 1141 (see figure caption 
for further references). It should be emphasized that all the 
determinations of stellar mass densities at redhsifts z > 2 
previous to this study have been made on pencil-beam sur- 
veys, with coverage areas < 5arcmin^. 

The analysis of the distribution of the stellar mass 
budget in Figure 1131 also allows us to extract an interest- 
ing conclusion: a minimum of ~ 45% and a maximum of 
~ 70% of the stellar mass at redshifts Zphot ^ 3 is con- 
tained in galaxies with assembled mass M > 1 X 10"Mq. 
To compute the lower limit (~ 45%), we considered the frac- 
tion of the stellar mass density assembled in galaxies with 
M > 1 X 10"'^"'^ Mq at {zphot) = 3.5 over the total stellar mass 



density at (zphot) = 2.5. The latter can be considered as 
an upper limit to the total stellar mass density assembled 
at (zphot) = 3.5. The upper limit to the fraction of mass 
locked in massive systems (~ 70%) has been obtained by 
directly comparing the stellar mass density corresponding 
to galaxies with M > 1 x 10^^ M© and the total value, both 
at (zphot) ~ 3.5. The large percentage of stellar mass locked 
in massive systems is related to the significant number den- 
sity of massive systems we find to be present before redshift 
Zphot ~ 4, and indicates that galaxy/star formation at very 
high redshift is an extremely efficient process. These mas- 
sive galaxies present at very high redshift have presumably 
been formed from the early collapse of the high peaks of the 
density fiuctuation field (e.g. Mo, Jing & White 1997). Such 
regions of high density are predicted to be strongly clustered, 
a fact which is perfectly consistent with the strong cluster- 
ing observed in red Kg-selected galaxies at 2 < Zphot < 4 
(Daddi et al. 2003) and in the ERGs up to the faintest mag- 
nitudes (R03). Evidence for a similar 'anti-hierarchical' be- 
haviour in favour of biased galaxy formation has also been 
found by other authors from the study of galaxies at redshift 
z ~ 1 (Kodama et al. 2004). Our results are consistent with 
a scenario, already suggested by previous studies, in which 
a significant fraction of the massive galaxies we see today 
have been assembled at very high redshifts (2 > 4), evolved 
to (I — Ks) ERGs by redshift « ~ 1 and possibly became 
members of massive galaxy clusters in the local Universe. 
This first detailed study of a significant deep sample of Ks- 
selected galaxies allows us to show and summarise in one 
plot (Figure IT^ ) the evolutionary line traced by the mas- 
sive Ks-selected galaxies and the massive ERGs since very 
high redshifts. 



5 PREDICTIONS OF SPITZER/IRAC 
MAGNITUDES 

HYPERZ provides the possibility of extrapolating the 
galaxy magnitudes beyond the observed bands using the 
best-fit SED template for each object. We used the observer- 
frame best-fit templates of each of our Ks-selected galax- 
ies to predict their expected magnitudes in three of the 
Spitzer/IRAC channels, i.e. 3.6/im, 4.5/im and 5.8/^m, as 
well as the derived IR colour distribution as a function of 
photometric redshifts. We present in this section the pre- 
dicted colour distributions for all our Ks-selected galaxies. 
The complete catalog of Ks-selected galaxies with the indi- 
vidual Spitzer/IRAC predicted magnitudes will be released 
in a future paper. To convert from fiuxes to magnitudes, 
we used the IRAC zero-fiux points: 277. 5Jy, 179. 5Jy and 
116. 6Jy at 3.6/im, 4.5/im and 5.8/im, respectively (Fazio et 
al. 2004). 

Figure \Wl shows the predicted Ks — [3.6pim] (a), Ks — 
[4.5/im] (b) and Kg — [5.8/im] (c) colours of Ks-selected 
sources versus estimated redshift. The filled circles corre- 
spond to the ERGs and the open circles to all the other 
Ks-selected galaxies. We observe that (I775 — Ks) selected 
ERGs do not have discriminating near-IR colours, and a 
pure IR colour cut does not seem to select the same popula- 
tion as an (I775 — Ks) cut, as was pointed out by Wilson et al. 
(2004) (these authors compared a pure IR colour cut with 
an (R — Ks) cut). From the analysis of Spitzer/IRAC data. 
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Figure 15. Predicted IR colours of Ks<22-selected sources as 
a function of estimated redshift. [3.6], [4.5] and [5.8] indicate 
Spitzer/IRAC magnitudes at the channels 3.6/im, 4.5/jm and 
5.8/im, respectively. 



these authors also conclude that a Kg — [3.6/im] red colour 
could be more effective to select high redshift (z > 1.3) 
galaxies. The predicted IR colours for our deep sample of 
Ks-selected galaxies leads us to a similar conclusion. How- 
ever, a Ks — [3.6/im] colour cut sufficiently red to avoid z < 1 
contaminants might miss some of the higher redshift objects. 
On the contrary, the Kg — [5.8/im] colour appears as a rather 
better redshift indicator, and a red Kg — [5.8/im] cut could 
more effectively separate the higher redshift galaxies. The 
relatively tight relation existing between the Kg — [5.8/im] 
colour and redshift suggests that the three very blue sources 
with Zphot ~ 2.5 — 3.0 might have erroneous estimated red- 
shifts. 

It is interesting to compare the predicted magnitudes 
for a subset of our sources with recent results obtained by 
direct photometric measurements on IRAC images. Yan et 
al. (2004) recently selected a sample of 17 IRAC objects with 
/„(3.6/im)//„(z85o) > 20 in the HUDF. We find that 8/17 
of these objects are members of our Ks<22-selected sample. 
The remaining objects are not in our sample either because 
they are fainter Kg-band galaxies or because they lie outside 
of our field of view^. 1/8 of the common objects (identifi- 
cation number 1 in Yan et al. 2004) has a very uncertain 
estimated redshift in Yan et al. (2004) and is excluded from 
their refined sample. The remaining 7/8 common objects 
are listed in Tabled Our estimated redshifts for the seven 
common sources are in good agreement with the Yan et al. 
estimations (median of Sz/{1 + z) Ri 0.07). We note that, 
however, our predictions underestimate the observed mag- 
nitudes in the IRAC bands by a median of ~ 0.5 mag. The 
comparison of the measured Kg magnitudes in both cases 
shows that this systematic offset cannot be attributed to, for 
example, aperture effects. However, given that the sources 
listed in Tabled have been deliberately selected by Yan et 
al. on the basis of their flux excess on the IRAC bands, it 
is unsurprising (and reassuring) that, for this highly biased 
subset of sources, our predicted IRAC fluxes are systemat- 
ically low. We anticipate a much better agreement of our 
predicted values with the observed IRAC magnitudes for 
the majority of our sample, when the GOODS/CDFS IRAC 
data become available. 



6 SUMMARY AND CONCLUSIONS 

In this work we have presented the redshift distribution and 
resulting cosmological implications for a sample of 1663 Ks- 
selected galaxies with Kg < 22 (Vega), and made an exten- 
sive comparison of the results with those obtained in C04 for 
the R03 ERG sample in the same field. This is the deepest 
significant study to date of Ks-selected galaxies and the role 
of ERGs within this population. 

We have studied the evolution of the rest-frame Ks-band 
LF and concluded that there is no evidence for a decrease of 
its bright end up to at least redshift (zphot) = 2.5. A similar 
conclusion has been reached for the ERGs in C04, because 
these objects account for most of the bright end of the Ks- 
band LF at high redshifts. Also, from the comparison of both 

^ The ISAAC field studied in this work only partially overlaps 
the HUDF. 
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Table 1. Ks<22-selected sources in the present sample which arc among Yan et al. (2004) IRAC extremely red objects. The columns 
list the following: 1) our identification number; 2) Yan ct al. identification; 3)-4) RA and DEC on the Ks-band images; 5) our estimated 
redshift; 6) Yan et al. estimated redshift; 7) our measured Ks-band magnitude; 8) Yan et al. measured Kg-band magnitude; 9)-10)-ll) 
our predicted magnitudes in the IRAC channels 3.6/im, A.bfim and 5.8/im, respectively ; 12)-13)-14) magnitudes measured by Yan et al. 
in the IRAC channels 3.6/ir7i, 4.5/i?7i and 5.8/im, respectively. All magnitudes are in the Vega system. To convert Yan et al. magnitudes 
quoted in the AB system to Vega magnitudes, we used: Ks^^^^ = Ks^g — 1.84, [3.6]vega = [3.6]ab ~ 2.79, [4.5]vega = [4.5]ab — 3.26 and 
[5.8]vega = [5.8]ab — 3.73. 
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the total Ks-band and the ERG LPs, we find that the latter 
appears to have a bell- like shape at redshift ZjjUot ~ 1, with 
little contribution to the faint end of the former. This result 
is consistent with other previous works, which determined 
that red and blue galaxies dominate different parts of the 
total LF in the local Universe as well as higher redshifts. The 
total Ks-band and the ERG LFs become progressively more 
similar with redshift, reaching a maximum at {zphot} = 2.0. 
The limiting magnitude of our sample does not allow us to 
properly explore the faint end of the LF above this redshift. 

We also study the evolution of the massive systems 
present in our sample of Ks-selected galaxies. We determined 
that a significant fraction (up to 20%-25%) of the galax- 
ies with mass M > 1 X 10" Mo is in place before redshift 
Zphot ~ 4. However, the assembly of most of the remaining 
massive systems appears to have occurred at later epochs, 
at redshift Zphot < 1.5. ERGs at redshifts Zphot ~ 1 — 2 ac- 
count for most of the massive galaxies which have plausibly 
evolved under passive evolution since very high redshifts, 
but 30%-50% of the galaxies with assembled stellar mass 
M > 1 X IO^Mq at redshift Zphot ~ 1 - 2 are not ERGs. On 
the other hand, ERGs account for only a fraction of the mas- 
sive systems present in the local Universe, indicating that, 
in any case, passive evolution could not entirely explain the 
construction of all the massive galaxies. 

We found that the comoving number densities of gala:x- 
ies with stellar mass M > 1 x lO^^M© at redshifts Zphot ~ 
1 — 3 is significantly smaller than the corresponding den- 
sities of dark matter haloes massive enough to host these 
systems. We suggest that this could be evidence for differ- 
ential timescales in halo merging and galaxy collapse. At 
redshifts Zphot > 3, however, the number of massive galaxies 
coincides with the corresponding number of haloes, proba- 
bly indicating that galaxy collapse was a much faster and 
efficient process at very high redshift, virtually simultaneous 
with the virialisation of the haloes. 

From the analysis of the distribution of the stellar 
mass budget at different redshifts, we conclude that be- 
tween 45% and 70% of the stellar mass assembled at red- 
shifts 3 < Zphot < 4 is contained in galaxies with mass 
M > 1 X IO^^Mq. These massive galaxies present at very 
high redshift have presumably been formed from the early 
collapse of the high peaks of the density fluctuation field. 



Our study of the evolution of massive systems up to red- 
shifts Zphot ~ 4 seems to confirm the existence of an evo- 
lutionary sequence in which most of the massive galaxies 
formed at very high redshift become (I — Kg) selected ERGs 
by redshift Zphot 1, and are the progenitors of the massive 
elliptical galaxies observed in local clusters. 

The follow up in the Spitzer/IRAC cliarmels will be nec- 
essary to directly observe the rest-frame Kg-band light of the 
high-redshift (z > 1) galaxies in our sample. This will allow 
an independent test of our extrapolated stellar luminosities 
and masses at high redshift. Moreover, the study of comple- 
mentary properties for Ks-selected galaxies (e.g. morphol- 
ogy) should shed further light on the nature and evolution 
of these objects. 
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